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A 6-ch()ice simple cubic lattice model of adsorption of an isolated polymer chain at a solution 
surface is investigated. The mean square components {x'^{N)) and (zmV)) of the end-to-end dis- 
tance are computed as a function of the adsorption energy per monomer unit in the limit of a very 
long polymer chain. In the calculation, one end of the polymer chain consisting of A^ monomer units 
is ccmstrained to lie in the surface; and (.r^(AO> and (zHN)) are, respectively, the mean square dis- 
placement of the free end of the chain parallel to the solution surface in one of the lattice directions 
and normal to the solution surface. The limiting value of {x-(N))IN as A^-^oc is a continuous func- 
tion of 0, the dimensionless adsorption energy per monomer unit, and is equal to 1/3 for 6 ^\n (6/5) 
and (1/2) [l+(l/4) (e*^ - \)-']-"^ for 6» ^ In (6/5). The hmiting value of {zHN))!!^ as A'-^oc is a dis- 
continuous function of 6 and is equal to 2/3 for 6^ < In (6/5), 1/3 for 6^ In (6/5), and for > In (6/5). 
The relation of these results to earlier investigations and the generalization of these results to other 
cubic lattice models is discussed. 

Key Words: Adsorption, chain polymer, critical energy, generating function, lattice model, parti- 
tion function, random walk. 

1. Introduction 

In this paper we extend the investigation of a random walk lattice model of polymer chain 
adsorption at a solution surface [1, 2| ' by computing the components of the mean square end-to-end 
distance in an adsorbed polymer chain. The random walk lattice model used is a generalization 
of the models introduced by Silberberg [3] and DiMarzio and McCrackin [4]. The physical proper- 
ties of the model reflect the competition between (1) the energy gained by the polymer chain occupy- 
ing as many surface sites as possible, and (2) the configurational entropy lost by the polymer chain 
remaining near the confining solution surface. In I and 11 the average fraction of monomer units 
which lie in the solution surface p{d, N) was computed as a function of the reduced energy of 
adsorption per monomer unit 6 = elkT. In the calculation one end of the polymer chain was 
constrained to lie in the surface, and the self-excluded volume of the polymer chain was neglected. 
In the limit in which the number of monomer units A^ is large, it was determined that there is a 
critical value of the reduced energy of adsorption per monomer unit 6c such that for > Or, the 
molecule exists in an adsorbed state. For example, lim ^{6, N) , the limiting average fraction 
of monomer units lying in the surface, is a positive constant independent of N indicating that a 
finite fraction of all monomer units lie in the solution surface. For d< 6c this limiting fraction 
is zero. 

The purpose of this paper is to calculate {x\N)) , (/(AO), and {z\N)), the components of the 
mean square displacement of the free end of the polymer chain as a function of ^ for A^ > 1 for the 
6-choice simple cubic lattice model. It was shown in I that the mean distance of the free end of 
the chain from the solution surface {z{N)) has the following discontinuous form for the 6-choice 

' Figures in brackets indicate the literature references at the end of this paper. 
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simple cubic lattice model 





{7Tibyi\ 


« 6* < In (6/5) 


lim (2(A0)/A'"-= { 


{2l3Tryi\ 


e = ln(6/5) 




0, 


e>\n (6/5). 



It is shown in section 4 that there is a corresponding discontinuity in lim (z^(N))IN whereas for 
the other components, lim {xHIV))IN = l\m. {y^{N))IN is a continuous function of 6. The recur- 
rence equations for the probability distribution of the free end of the polymer chain are formulated 
in section 2 and are solved formally in section 3. Explicit expressions for the second moments 
are evaluated in the limit A^^ ^ in section 4. The results are summarized in section 5 and the 
qualitative behavior for the 2-choice and 4-choice simple cubic models is deduced. 

2. Recurrence Equations for the Probability Distribution of the 
Free End of the Polymer Chain 

We consider a simple-cubic lattice model of the solution-surface system in which the solution 
surface corresponds to the x-y lattice plane through the point z = i). Successive lattice planes 
through z= 1, 2, . . . represent the solution phase of the system. Polymer chain configurations 
in the solution correspond to paths generated in a random walk on the lattice between nearest- 
neighbor sites. The physical presence of the surface is introduced by considering only random 
walks in the lattice planes through z = (). 1,2, . . . which never enter the lattice plane through 
z = —l. In the absence of a solution surface, all random walk configurations of a given length are 
equally likely. However, we are primarily interested in the influence of an adsorbing solution sur- 
face on the average conformation of a polymer chain. In our lattice model, all random walk paths of 
A^ steps with n steps lying in the surface layer z=i) have the same a priori probability. Relative 
to a random walk configuration of A^ steps with n — 1 steps lying in the surface, the a priori proba- 
bility of a walk with n steps in the surface layer is greater by the factor e^ where = elkT and e is 
the adsorption energy of a monomer unit. 

For convenience, we use a random walk terminology. Consider the problem of computing, 
for a random walk originating in the surface layer, the unnormalized or relative probability 
P{x, y, z; A^-h 1) that at the A^+ 1th step the random walker is located at lattice point (x, y,z), where 
z ^ 0. The relative probability P{x. y, z;A^+r) is related to the relative probabilities at the 
A^th step by the relations 

1 , (1) 

P{x.y.z'.N+\)=-{E-,^E--^E^,^E-y+Et+E-,}P{x,y,z;N), z^l 

and PU, y, 0; N +\) = ^e'{Ei + E-^ E^^-^ E-^Et]P{x, y, 0; AO (2) 

where E^ are operators defined by the relation 

E^P{x, y, z; AO = PU ± 1, y, z; AO. (3) 

The operators £'J and E~ have similar definitions. Equation (I) describes the relation between 
the relative probabihty of being at lattice site x, y, z with z ^ 1 at the A^-h 1th step and the relative 
probabihties of being at neighboring sites at the Mh step. The factor e^ in eq (2) accounts for 
the fact that relative to those configurations where z ^ 1 at the A^+ 1th step, the relative probabil- 
ities for those configurations where z = at the A^-hlth step are greater by the factor e^. The 
absence of Ez in eq (2) is related to the fact that the random walker enters the z = layer only 
from one direction. Equations (1) and (2) will be solved in section 3 for the initial condition 
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P(x, r, z:0)= (4) 

[0, for all other lattice points. 

From the solution, the mean scjuare end-to-end distance after A^ steps can be computed from the 
expression 

K?=MAO) + (r(/V)> + (z-(A/))= I; f^ ^(x' + f+z')P{x,y,z-N)jQ(e,N) (5) 



''^'""' Q{e,m= X 2 ^Pix,y,z;N). (6) 

3. Solution of Recurrence Equations 

The recurrence eqs (1) and (2) can be solved for the starting condition eq (4) by introducing 
a generating function as in 1 and U. It is convenient first to rewrite eq (2) as 

[1 -a-e')]nx. y. 0; A^+ 1) = (l/6)[^^ -f £; + ^; + ^^ +£^]P(x, y, 0; M). (7) 

Multiply eqs (1) and (7) for P{x, y, z; /V+ 1) by (27t)~'^'~ exp {ix^-\- iyr] + /z^} and sum over all integer 
values of x and y and all z ^ 0. The result is 

-(1-c ^)(27r)^^/^/;(^, 7]; 0; /V + l) + r;(^, t], ^; /V+1) 

-^ [cos f + COS 7) + cos ^]r;(^, r;, ^; AO"^ c-'^(27r)-'/^p(^, r?; 0; AO (8) 



X X X 



where 6'(f , t], ^; AO = (277)--y2 ^ S S ^'''' >■' ^' ^ ^'^P Hx^+iyV + f^D (9) 

J'=-x //=-x 2 = 



x X 



and i>{i,'n;z; IV) = (271)-' ^ 2 ^(-«. r, 2; AO exp (i%^+ jyrj). (10) 

X=-oc //=-x 

Next multiply eq (8) by u/^^^ and sum over all values of A^ to obtain 
n^, 7), ^; w)-G{^, J], e; 0)-(l-e-^)(27r)->/^[p(^, t?; 0; m;)-/X^, t?; 0; 0)] 

= -z^[cos ^ + cos T7 + COS £]r(f, 17, ^; z^)— -e-'K277) ^/-p(^, 17; 0; w) (11) 
where Hf, 7?, ^; m;)= |] 2i;'r;(^, 7], ^; AO (12) 



and 

A'=0 
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p(^, 7]; z\w)=^ w^p{^, 7]; z\ N). (13) 



In eq (11) the quantities G(^, 17, (,: 0) and p(^, 17; 0; 0) are simply 

G(^, T^, ^: 0) = (27r)-^/V 



(14) 



and 



p(^, T?; 0; 0) = (27r)-ie^. 



(15) 



Substituting eq (14) and (15) in eq (11) and solving for H^, 17, ^; w), one obtains 



(27r)-'^/2 + (2^)-i/2 



n^, T], ^; zi')^ 



1 — e ^ — - we '^ 
o 



p(^. 17; 0; w) 



1— -z^'[cos ^ + cos 17 + cos ^] 



(16) 



E([uation (16) is an implicit ecjuation for p(^, 17; 0; w) because according to the definitions of 
n^, 7], ^; w). G(^, J], i\ /V), p(^, r?; z; w;), and p(^, t?; z; /V) in eqs (9), (10), (12). and (13) 



X XXX 



i^(^, 1). C: iv)={27t)--^I'' 2 2X2 ^'^-^- >'• ^^ ^^^"^^ ^^1^ ( Jx^+ /yr?+ /2^) 

.r=-x (y=-x ^=0 A'=() 



= (277)-^/^ ^e'^^pii, 7]\ z; 2^). 



(17) 



The following relation for determining p(f , 17; 0; z^) can be obtained by multiplying eq (16) by {'Irr) ^^~ 
and integrating with respect to ^ from — tt to tt 



p(^, T]: 0: w) = {27t)-'U^. t); m;) + 



[1 -c-%)(f, rj: w)--^ivl-^^. 17: zi') p(^, 17: 0; zi;) (18) 



where 



liA^. T?: z^) = (27r)- 



e^'^y^ 



1 — z^(cos ^ + cos 17 + cos 0/3 



1 — z^(cos ^-\- COS T7)/3 



wl?> 



1 — z^(cos ^ + COS 17)/ 3 



-1/2 



1- 



w/S 



1— ^;(cos ^ + cos T7)/3, 



1/2 1 |A-| 



wis 



1— m;(cos ^-hcos T7)/3 



(19) 
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Solving eq (18) for p(f , rj: 0; w) and substituting the result in eq (16), we obtain an explicit solution 
of the recurrence eqs (1) and (2) for the initial condition (4) 



ri^,r^,C;w) 



1-3/2 



I l+(W6)[/,(g, 7): w)-e-^^Io{^, r]; w)] \ 

1 -Mcos ^ + cos T7 + COS ^)/3 l+(/tV6)/,(^,r;; w)-{\ - e-^)Io{^, r/; w)\ 



(277 



(20) 



The mean square end-to-end distance (5) can be expressed in terms of the generating function 
n^, T/, ^; w). First note that the normalizing sum Q{6, N), defined in eq (6), is proportional to 
the coefficient of u/"' in the expansion of HO, 0, 0; w) [see eq (17)] 



X X 



r((). 0. 0; w) = {27t)"'"'' X ^'' 2 S S^'*-^' y- ^^ N) = i^2^)-'"-^ ic'QiO. A) 



(21) 



Usin^ (>auc*hy's formula and ecj (20), we obtain the following explicit expression lor (J{6, N) 



277/ Jr-, W^^^ 27TI JCu 



l + [(l-?W3)/(l-z^)]^ 



dw 



[{l-w){l-wl3)y'''-\-{2wl3)^2e-' lu^ 



(22) 



where 60 is a counterclockwise contour of integration around w = {). It aisc) follows from v(\ (17) 
and the symmetry in the x and v coofdiiiatcs thai ihc numcraloi" in v(\ (S) is propoit ional lo the 
coefficient of iir' in 






m. 0, ^; Hi 



■X. XXX 



= (277)--'/2 2 «^'' 2 2 2 <^' + >'' + ^'' ''•^- ■>'• ^= ^'- (23) 



The explicit expression for the numerator in eq (5) is 

-(277)"'^/2 



2JAN)^^z{N)-- 



277/ 



'i\' 



2ii + li 



r(f, 0, ^; m;) 



dw 



(24) 



^ = 0,^ = 



or 



./r(/V) = 



(/^/3)^ 



1 m;/3 1 



277/ [(1-m;/3)i/2( 1-^^)3/2 D((9, m;) 1 - u; Z)'^(0, w;) 



1 + 



1 - W \l/2 



1-m;/3 



[1 - (2W3) + (1 - W3)i/2(1 - ^)i/2] j-jl^ (25) 



0im=^\ _iZ3_(ll^W3)^W 

'^ ^ 27TiJc,D{e, w) \ {l-w)''l^ Jw'' 



(26) 
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"^^^^^ D{e,w) = [{l-w)a-wl3)yi^-l-i2wl3) + 2e-'. (27) 



In eq (24) J'j>{N) and J'ziN) are, respectively, the unnormalized mean square % and z components 
of the position of the free end of the chain. The exphcit expressions (25) and (26) for J^j.(A^) and 
J^zi^) are evaluated in the limit A^^ ^ in the next section. 



4. Evaluation of /? v for iV^ 1 

We now consider the problem of obtaining asymptotic formulas in the limit A^ ^ 1 for the 
X and z components of R\ in eq (5). The explicit expressions which must be evaluated are 

and 

where ^AN), J^ziN), and Q{e, N) are given in eqs (25), (26), and (22). The details of the evalua- 
tion of these contour integrals in the limit A^ > 1 are straightforward but lengthy and depend as in 
I upon whether 6>6c, 6 = 6c, or < Oc where ^c = ln (6/5). All three integrands have branch 
points at w=l and w = 3. As in I, we introduce a cut in the w-p\ane between these branch points 
as shown in figure 1. In addition to the branch points, when ^> 6> 0^ there is a zero of D{0, w) 
on the real axis between w = and w=l. Details are given in I where it is shown that the expres- 
sion for this zero is 

u;+ = 6{-2(l-e-^) + [4(l-e-^)2 + e-^(l-e-^)]i/2}. (28) 

It is also shown that w+ approaches as ^ ^ oo and it approaches 1 as ^ ^ ^c = In (6/5). Thus when 
0= Oc, the pole of the integrands at w+ coincides with the branch point at w=l. In the remainder 
of this section we obtain asymptotic formulas for the contour integrals (25), (26), and (22) in the 
hmit A^ > 1 for the three cases 0> Oc, 0=6c, and < Oc. 

9 > 9r- For e>ec the contour integrals J^AN), ^z{N), and Q{e, N) all satisfy the relation 



^/c/^-^^-^-2^L^^^-^^-^2^L/(-)^- 



(29) 



where the integral on the left-hand side of eq (29) stands for J' xW, ^ z{N) , or Q{0, N) and the 
contours are shown in figure 1. In the Hmit N>1, the second integral on the right-hand side of 
eq (29) is negligible compared to the first integral. Therefore, the asymptotic values of -J^j^(7V), 
J" z{N), and Q{d, N) are given by the residues of the appropriate integrand at w+. The values are 



Q{e,N)=-w^-^''^'^{\^[[i-wA^)l{i-w^)yi^}ID'{e,w^) oo) 



^z{N)=-\w^-^\\-w^lzy^^l{\-w^Ym'{e,iv^) od 
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Figure 1. Integration contours Co, Ci, and C2 in the complex w-plane. 

The contour C-i is indented around the left end of the cut terminating at w= 1. 



'AN) 



N iv+13 



1 + 



I— W i 

1 — UJ+I'S 



1/2 



x[l-{2wJ3)^{l-wJ3)'IHl-w^)'l'] [D'ie, w^)y' (32) 



where D'{6, w+) denotes {dldw)D{6, iv)\ir=ir and where terms proportional to A^^ in eq (32) have 
been ne^ilected. Equation (30) is a speeial case of eq (37) in I. Combining eqs (30) and (31) 
we obtain for the mean scjuare z-coniponent of (Ry) 



(z'iN)} 



^ 3 il-w+fl'' 



1 + 



1 — w-^13 

1 —w+ 



1/2 






1 — W7+/3 



ii\ 



1/2 



(33) 



Thus, it is seen that if 6 > hi (6/5), then iv+ < 1; and the Hmiting value of (z-^{N))IN is 



hm {zHN))IN=0. 



(34) 



Conil)inin«: etjs (30) and (32), we obtain for the mean square x-component of {R\) 



{xHN)} 



N 



+ 



1 - (2m;+/3) 
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[D'ie^w^)]-'^- 



1 e-^ 



4 1 - e-^j 



-1/2 



(35) 



In obtaining eq (35) we used the relation 



D'iO. w^] 



-2-^ w+ 



l{\-w^yi^' {i-w^ld>yi' 



Ll+(2W3)-2e ^ 



(36) 



and eci (28) lor w^. The hmiting value of {x\N))IN is 



lim (^.(AO)/yV=-l,+-^j 



(37) 



5=^c. For 6=ec=\n (6/5) the factor [D{e, w)]"' appearing in the intejirands oiQ{0, N) . 
Jx{N). and JAN) can be written as 



[D{(),..w)\-' = 



S + w 



l-w/SV^- 2 



(38) 



When [D(dc, w)]"' is replaced by eq (38), the expressions for QiO, N), Jz(N), and .^xiN) become 



Qie.N) 



[I Jc o 



'Ztti J c 5 + IV 



l-w/3 , 1 [l-w/3 



] — w 3 \ 1 — it' 



"^ 2] dw 

3 J m;-^-^'' 



(39) 



J^.(AO = 



1 

"2m 



/ 



5 + if 



1-m;/3 2 (1-W3)'/^ 
(l-wF 3 (1-m;)'"-^ J 



dw 



(40) 



^x(/V) = 



277^ J^jl S-^W 



1 



2/3 



(l-Z^)i^ (1-^3)1/2(1-2^) 



3/2 



+ 



27z. 



(5-f w;)- 



1-W3 4(1-W3) 



1/2 



(1— m;)- 3 (1— z^ 



,3/2 



- + 



4/9 1 
l-itJ 



1 + 



1 



1-W.V3 



1/2 



[1 - (2m;/3) + (1 - u;/3)i/Ml-w;)'^'] 



G^zi; 



(41) 



The three contour integrals now satisfy the relation 

- — : R{w)dw = - — : R{w)dw. 

2771 Jc\. 27TI Jc\, 



(42) 
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It is readily verified that in each case the dominant contribution to the integrals (40) to (42) around 
the contour C2 comes from the most singular terms in (1 — z^). In (40) to (42) these singularities 
are simple poles. Evaluating the residues at these poles and retaining terms proportional to the 
highest power of A^ only, one obtains 

i){0,:N)=\. (43) 

y,(/V)sA'/3. (44) 

and J'AN)=Nli. (45) 

The limiting values of (zHN)}IN and {xHN))IN are 

lim (f{N))IN = \ (46) 

and lim {xHN))IN = \' (47) 

$<$<• For 0<0r. the ccmtour integrals (22), (26), and (25) satisfy the relation (42). The 
dominant contribution to the contour integrals for Q{0. N), -/^(A^), and -//(A^) around Cj comes from 
the most singulai- tcinis in (1 — w) in the imnic(hatc vicinity of w=}. The expressiims lor these 
dominant ('()ntril)uti()ns are 



<— i/:(va""i^iST. 



Ztti 



<■" {l-wl3yi^ 1/3 dw 
a {l-wf- D(fi. If) ur'" 



(49) 



^"^ "^-^'^ ~ 2m J, (1 - wISy^Hl ~ wfl^ \D(e, w) ) J^> ^^^^ 



where the points Q and R are shown in figure 1. It is possible to avoid the cxj)ncit cvahiation oi" 
the contour integrals in (48) to (50) because we are interested only in obtaining asymptotic fornuilas 
for yAN)/Q{0. N) and yAN)IQ{0. N). To demonstrate this fact, integrate by parts in (50) to reduce 
the order of the singularity in (1 — w). The result for -/r(AO is 

^J_ P 2NI9 } dw 

'■* ' ~ 277/ Jg (1 - wyl^' (1 - w/SyiWid, w) mA ^^1) 

where we have neglected contributions of lower order in A^. We also remark that the singularities 
in l — w in (49) and (50) are identical, as are the singularities in (l — w) in (48) and (51). Conse- 
quently, the asymptotic formula for ^zW can be deduced from that of .yj^(AO by multiplying -/r(AO 
by the factor 

hm 3(1-m;/3) = 2; 
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and the asymptotic formula for J^xiN) can be deduced from Q{0, N) by multiplying Q(6, N) by the 
factor 



lim 



2N 1 N 



w^\ 9 l — wl3 3 



Thus we have the relations 



yAN)=NQ{d,N)l3 



(52) 



for N>1. The hmiting values of {xHN))IN and {zHlV)}IN are then 



and 



hm {xHN))IN=ll3 
lim {zHN))IN=2l3. 



(53) 
(54) 



5. Summary and Remarks 

The limiting mean square components of the end-to-end displacement parallel to the solution 
surface are equal and have the value 



Hm {xHN))IN= Vim {f(N))IN 



3' 



3^ 

I 

2 



^ ^ < In (6/5) 



^ = ln (6/5) 



l-f7(e^-l)-^ 
4 



(55) 



6^ > In (6/5). 



Thus lim {x\N))/N is a continuous function of the reduced adsorption energy per monomer 
unit which increases monotonically from Va to V2 in the ^-interval [In (6/5), 00). The hmiting mean 
square component of the end-to-end displacement normal to the solution surface is 



hm iz\N))IN-- 



2 
3' 

1 
3^ 



0^6/<ln(6/5) 
^ = In (6/5) 

^ > In (6/5). 



(56) 



It is seen in (55) and (56) that at 6=0c the limiting mean square components of the end-to-end 
distance are equal to 1/3, the value which would be obtained in the absence of the solution surface. 
This behavior at the critical energy was deduced by DiMarzio and McCrackin [4]. DiMarzio [5] 
has obtained our result in the special case ^=0, but for a general lattice model. In this case 
(6^ = 0), the mean square normal component is twice the mean square lateral components which 
are in turn equal to the value obtained in the absence of the solution surface. It is seen in (55) 
and (56) that as 6^-^ ^ and the molecule becomes confined to the solution surface, the mean square 
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Tablk 1 



Values of mean square components of the end-to-end displace- 
ment in the 6-choice, 4-choice, and 2-choice simple cubic 
lattice models 



e=er 



{xHN))IN 



1/3 
1/3 

1/2 



{z'(N))IN 



2/3 
1/3 




Figure 2. Limiting form of the mean square end-to-end 
distance per step versus the adsorption energy 6. 



lateral components approach the value which would l)e obtained for a random walk on a 4-choice 
[)lanar square lattice, (^iomhinin^ the limiting forms for the mean square components of the end-to- 
end distance, we obtain finally for the mean square end-to-end distance per step 



im R'^IN-- 



1, 



^ ^ < In (6/5) 



^ = ln (6/5) 



I 



l+4(e^-l)- 



(57) 



<9>In (6/5). 



This function is plotted versus in figure 2. 

We now consider the question of determining qualitatively the components of the mean square 
displacement in the 4-choice and the 2-choice simple cubic lattice models. These models have 
been studied [1, 2] and the critical energies 6c = \n (5^^^ — 1) and (1/3) In 2, respectively, have been 
obtained. We expect that analogous to the result in the 6-choice simple cubic lattice model, the 
mean square components of the end-to-end distance in the adsorption energy range < ^ < ^r are 
equal to the values obtained at ^ = by DiMarzio [5]. We also expect that at 6= 6c, the mean 
square components of the end-to-end distance are equal to their free space values. As 0^^^ 
the mean square components parallel to the surface approach the values for the corresponding 
planar lattice. In table 1 we list the limiting values of the mean square components of the end-to- 
end distance which are obtained on the basis of the foregoing arguments and note that values listed 
are the same in the 2-, 4-, and 6-choice simple cubic lattice models. Only the critical energies 
and the rates of approach to the 6 = ^ values are different. The results for the 4-choice simple 
cubic lattice model are useful as reference values in McCrackin's Monte Carlo calculations of 
/?|r taking into account the self-excluded volume of the polymer chain [6]. 
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